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A  I 

For  a  randomized  complete  block  design  with  additive  block  effects, 
an  asymptotically  distribution-free  simultaneous  confidence  region  of 
pairwise  treatment  differences  is  presented.  The  corresponding  confidence 
bound  has  an  explicit  form  and  is  easily  obtained.  An  example  is  provided 
for  illustration  purpose.  The  case  of  treatment  against  control  is  also 
discussed. 
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1 .  INTRODUCTION 


Suppose  that  K  >  2  treatments  are  applied  once  each  to  n  different  blocks. 
Let  be  the  response  of  the  ith  treatment  in  the  jth  block 

(i=l,  j=l,...,n).  The  model  often  used  for  this  experimental  setting  is 

the  linear  model  in  which  the  observations  X, .  can  be  written  as 

ij 

(1.1)  Xij  *  y  +  a±  +  +  e^. 


K 

where  the  a's  are  the  parameters  of  interest  (treatment  effect),  E  o  =  0, 

i=l  1 

tf's  are  nuisance  parameters  (block  effects)  and  e_.  =  (e^  , . . .  )  ' ,  j=l,...,n 

are  independent  and  identically  distributed  random  vectors  having  a  continuous 
joint  distribution  function  which  is  symmetric  in  its  K  arguments  (  This  relaxes 
the  conventional  assumption  of  having  independence  and  identity  of  distributions 
of  all  the  nK  error  terms.). 

Oftentimes  a  global  test  for  =  . . .  *  oR  is  of  less  interest  and 

one  may  feel  that  the  simultaneous  inference  on  the  pairwise  treatment  differences 
-  a^,,  l£i<i'<K  is  more  desirable  (c . f .  Miller  (1966,  1977)).  There  are  several 
nonparametric  pairwise  multiple  comparisons  procedures  available  for  this  case 
(c.f.  Puri  and  Sen  (1971),  Hollander  and  Wolfe  (1973)  and  Hettmansperger  (1975)). 
However,  they  either  only  utilize  the  intrablock  comparisons  and  have  low 
efficiency  or  involve  complicated  inversion  procedures  for  obtaining  simultaneous 
confidence  bound  of  pairwise  treatment  differences.  Although  Sen  (1969)  has  pro¬ 
vided  a  simultaneous  confidence  bound  to  -  a^,  based  on  two  sample  Chernoff- 
Savage  rank  order  statistics,  the  derivation  of  his  procedure  is  not  obvious 
(see  Puri  and  Sen  (1971),  p.  331),  and  the  bound  he  obtained  is  not  in  an  explicit 
form  so  that  numerical  method  is  sometimes  required. 

In  this  article,  we  utilize  the  information  contained  in  interblock  comparisons 
and  provide  an  asymptotically  distribution-free  simultaneous  confidence  region 


2 


of  pairwise  treatment  differences.  The  corresponding  bound  has  a  simple  and 
explicit  form  and  can  be  easily  obtained.  The  case  of  treatment  against  control 
is  also  discussed.  An  example  is  presented  for  illustration  purpose  in  Section  3. 


2.  ASYMPTOTICALLY  DISTRIBUTION-FREE  SIMULTANEOUS  CONFIDENCE 
REGION  OF  PATRWISE  TREATMENT  DIFFERENCES 

To  eliminate  the  nuisance  parameters  B's  in  (1.1),  we  consider  the  aligned 

observations  Y  =  X  -  X  where  X  is  a  symmetric  function  of  X.  X„.  , 

ij  ij  •  J  •  j  ij  KJ 

such  that  X  +  a  is  the  same  function  of  X,  +  a,...,X  +  a  for  all  -°°<  a  < 

•  J  lj  kj 

— 

Typical  X  .  are  the  block  average  (X  .),  the  median  of  X  , . . . ,X„  ,  the  Winsorized 
•J  -J  ij 

or  trimmed  mean,  etc.  In  this  article,  we  let  X  be  X  .  so  that  Model  (1.1) 

•  j  •  J 

can  be  rewritten  as  Y^  “  “i  +  eij  » •  •  •  >K;  j  =  1 . n)  .where  e  *  e^  “  e  j 

and  e  is  the  jth  block  average  of  e . e  .  It  follows  from  the  inter- 

•J  J-j  KJ 

changeability  of  e^,...,e].j  that  the  distribution  function  of  E^j » •  •  •  »eKj  » 
j=l,...,n  is  symmetric  in  its  K  arguments.  Let  the  marginal  distribution 
function  of  e^  be  G. 

Now,  define  a  scoring  function  <{>  for  comparing  Y^  and  Y^,^,  by 


(Yij’  = 


Y4J  >  Yitj, 


-1’  Yij<Yi*j' 


Let  p..,  =  E<J)(Y.  .,  Y. ,.,),  where  i^i'  and  j^j'.  Under  H  ,  p  ,=  0.  Consider 
n  ij  i  j  0  ii 

n  n 

U  =  Z  l  (4> (Y  ,  Y  ,  ,)  -  p,  ,),  1  S  i  <  i'  S  K,  which  is  the  usual 

11  *n  j  »=i  1J  1  J  11 

Mann-Whitney  two  sample  statistic  (though  based  on  matched  samples) . 

Several  lemmas  are  needed  to  derive  the  simultaneous  confidence  bound  to 

al-«i,,  H  i  <  i'  a. 
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Lemma  1.  As  the  random  vector  <n  U.t,  / 2  >  converges  in  distribution 

-  li  ,n 

to  a  K(K-l)/2  dimensional  normal  random  vector  with  mean  0  and  covariance  matrix  $. 


(Proof) . 

Lemma  2 . 


See  Appendix  A. 

n~  '2U 

Under  H-,  the  asymptotical  covariance  structure  of  < - r - : — - — > 

0  2(o2-t)/* 


identical  to  that  of  the  vector  having  K(K-l)/2  components  Z^-Z^',  l^i<i'sK, 
where  (Z^,...,^)'  ~  NK(0,1K),  1R  is  the  KXK  identity  matrix,  a2  =  E(G(Y^))2  =  Va 
and  t  =  E(G(Yij)G(Y1,j)). 


(Proof).  See  Appendix  A. 


Lemma  3.  Under  H^,  */6  £  T  5  V3  • 

(Proof).  Using  the  fact  that  )cov(G(Y.j,  ),  G  (Y^ ,  ^ ) )  |  ^  (var  G(Y^  )  var  G^.^))  2 
=  V12  »  the  proof  is  straightforward. 

We  remark  that  the  bounds  of  t  in  the  above  Lemma  are  attainable  (c.f. 
Hollander,  Pledger  and  Lin  (1974),  p.  180). 

Lemma  4.  Under  H  ,  lim  P(  |  U±i,  n  |  <  2q^  n  /3//6,  l<i<i’<K)  2:  1-y ,  where 

n*»  ,n 

0 <  y  <  1,  qy  is  the  lOO(l-y)  percentage  point  of  the  distribution  of  the  range  of 
K  independent  unit  normal  random  variables. 


(Proof).  It  follows  directly  from  the  above  lemmas. 


Lemma  5.  For  i<i’ ,  suppose  that  the  differences  ^  # ^ ,  are  distinct, 

ii*  if 

j , j '=1, . . . ,n.  If  D  <  ...  <  D  9  denote  the  ordered  differences  Y.,-Y.,.,, 

(1)  ( n )  J 

then 

iH  n  n  2 

D(l)  S  “i-0!'  lf  and  °nly  if  2  ^  (Yij _ai »  y i » j  » “a i » )  <  n  -21 
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and 


(m) 


ai~ai * 


if  and  only  if  £ 

j=l 


Z  ♦(Ylj-ai,  Y^.-c^,)  *  n2-2nrt-2. 


(Proof).  See  Appendix  A. 


Now,  we  are  ready  to  present  a  simultaneous  confidence  bound  to 
a^-a^, ,  l£i<i'<K. 

Theorem  1.  lim  P(D^.  £  o^-a^,  D(m)  *  ^  1~Y , 

n  -*»  '  *  '  ' 

where  if  *  [n  72  -  qj  n  //6],  m*[n2/2  n  ^  //6]  +  2,  and  [•]  is  the 

K  K 

greatest  integer  function. 


(Proof).  It  follows  from  the  above  lemmas. 


So  far  we  have  assumed  that  there  are  no  tied  observations.  When  there 
are  ties.  Lemma  5  will  no  longer  be  valid.  If  in  practice  the  ties  are  the 
result  of  rounding  to  the  nearest  multiple  of  e,  some  modifications  can  be  made 
to  guarantee  the  validity  of  Theorem  1.  Let  the  original  responses  giving  rise 

I 

to  the  (rounded)  observations  be  for  which  Model  (1.1)  is  appropriate, 

then  |X..  -  X.  |  <  %  and  hence  |X  -  X  .  |  5  E/? ,  where  X  is  the  jth 

-tj  lj  »J  *J  c  *J 

i  * 

block  average  of  X^  , . . . ,  Xj^ j  .  It  follows  that 


(2.1) 


'Yij  "  Yijl  S  e>  where  Yij  *  Xij  “  X.j* 


»  »  jlA  * 

If  the  ordered  differences  Y. ,  -  Y. are  denoted  by  E  ...  E  9  » 

3  3  (1)  (n  ) 

Lemma  5  holds  when  D  replaced  by  E.  However,  from  (2.1),  |  S  2  c. 

ii1  ii' 

Therefore,  if  II  [E,..,  E,  .)  is  the  1-y  simultaneous  confidence  region 
l<i<i'<K 


I 
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of  a  -a  , ,  l<i<i'<K  in  Theorem  1,  then  II 

1  1  l<i<i'<K 

a  1-y  simultaneous  confidence  region  of  a^-a^t • 


ii ' 

2c,  D,  .  +  2c) 
(m) 


is  also 


3 .  AN  EXAMPLE 


We  present  a  numerical  example  in  this  section  for  illustration  purpose. 

The  data  in  Appendix  B  were  obtained  by  Woodward  (1970)  to  compare  three  methods 
of  rounding  first  base  to  reach  the  second  base.  The  three  methods,  "round  out", 
"narrow  angle",  and  "wide  angle"  are  illustrated  in  Hollander  and  Wolfe  (1973, 
p.  142). 

Each  entry  in  Appendix  B  is  an  average  time  of  two  runs  from  a  point  on  the 
first  base  line  35  ft.  from  home  plate  to  a  point  15  ft.  short  of  second  base. 
Here,  players  are  blocks  and  methods  of  rounding  first  base  are  treatments. 

The  observations  were  rounded  to  the  nearest  multiple  of  e  =  .01.  For  error 
probability  y  =  .1,  we  obtain  i  =  119  and  m  ■  365  from  Theorem  1.  It  follows 
that  a  90%  simultaneous  confidence  region  of  a^-a^,,  lsi<i's3  is 


rD 


119  "  2e’  °365  +  2e)  *  tDU9 


which  is  L-.07 ,  .09)  x  [0,  .19)  x  [0,  .15). 


4.  REMARKS 

Suppose  that  Treatment  1  is  a  control  and  the  rest  (K-l)  treatments  are 
under  investigation  as  possible  improvements.  Then  by  the  same  argument  as  we 
gave  before  for  all  treatment  comparisons,  a  1-Y  simultaneous  confidence  region 
of  a^-cx^,  2<i^K  can  be  obtained  as  follows. 


Theorem  2.  lim  P(D 
n  ->■  °° 

2  Y 

!.  =  [n  /2-  £  n 
K-l 


il 


tt) 

3/2 


s  a. -a.  <  D,  .,  2siSK)  2  1-or,  where 
t  i  fra; 

//31,  m  =>  [n  2/z  +  n3/2//3]  +  2 

K-l 
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Y 

and  £  is  the  upper  Y  percentage  point  of  the  maximum  absolute  value  of 
K-l 

(K-i)  N(0,1)  random  variables  with  common  correlation  %  (c  f.  Hollander  and 
Wolfe  (1973),  Table  A.  14). 


APPENDIX  A 


,0 


Proof  of  Lemma  1.  Let  $  ,  (y)  =  E<f>(y,  Y.  ,  . ,)  -  p  .  ,  and 

ii  •' 


ii' 


(y)  =  K<f>  (Yij  ,y)  -  piif.  Also,  let 

g(ViJ,Yi..,)  -  ♦<Ylj>Y1.j,)  -  PU,  -  ♦?1.«iJ)  -  *it.  Pf1>v>  ---I 

Uli '  ,n  *  ^  <V  +  *ii'<Yi'J)>  '  2  J21(C<IlJ'“1'>"G<ll'j'“l)",>ll'>' 


Then 


_  3/  *  2  -3  n  n  n  n 

(A-l)  E(n  *  U,,t  n  -  n  2  U...  )  =n  X  El  E  h(j ,  j  ’  ,k,k' )  , 

11  *n  11  ,n  j=l  j'=l  k=l  k'=l 


where  h(j,j',k,k')  -  ELg(Yi^ ,Yi,j ,)g(Ylk,Yi,k,)  ].  Because  g  is  bounded,  we  can 

3 

ignore  any  u  terms  of  the  sum  in  (A-l)  if  u  is  of  order  o(n  ).  Consider  the 

following  cases  for  which  the  number  of  terms  is  with  order  larger  than  or  equal 
3 

to  0(n  )  (where  j,j',k  and  k'  represent  four  distinct  indices): 


(1)  h( j , j ' ,k,k' )  =  Eg(Ylj,Yi,jl)  Eg(Ylk>Yitk,)  =  0; 

(2)  h(j,j,k,k’)  h(j,j\k,k)  =  0; 

(3)  h(j ,j ' ,j ,k')  =  h(j,j',k,j')  =  E{E[g(Yi.jYi,jl)g(Yij>Yitk,)|Yij]} 

-  E{E[g(Ylj,Yi,jf)|Yij]  E  [g(Yij,YiIk,)|Yij]j  -  0; 

(4)  h(j  ,j  '  ,k,j)  -  h(j,j’,j',k’)  =  EE[g(Yi.,Y1,j1)g(Ylk,Yi,j  >1*^  *^2 

-  E{E[g(Ylj,Ylfjt)|Ylj,Y1,j]ECg(Ylk,Ylfj)|Ylj,YilJ3}  -  0. 
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-  3A>  -V„  *  0 

It  follows  that  E(n  U, . ,  -  n  *  U..,  )  ■*  0,  as  n-*-“>.  By  Corollary  6  of 

li  ,n  n  ,n 

m  T/  ^ 

Leluuann  (1975,  p.  289)  and  the  fact  that  the  vector  <n  *  U,Jt  /2  >  has  a 

li  ,n 

K(K-l)/2  dimensional  normal  limiting  distribution  N(0,$),  the  vector 

_  3/,, 

.  '  ‘  II  /  n  V  1- _  .  1  t  1  >  •  .  •  It.  ••  .  .  „  1  \  n  __ 


n/2  >  has  the  same  normal  limiting  distribution  N(0,$). 


Q.E.D. 


Proof  of  Lemma  2.  Under  the  H  :  a  =  a  =  ...  =  a  Y  ,  ....  Y  are  inter- 

U  X  ^  K.  X  j  K,  j 

changeable  and  have  a  common  continuous  marginal  distribution  function  G. 

A  typical  element  o^,^,  of  %  is  E^Y^  )-G(Y .  ,  J)  (G(Ykj)-G(YRl ,  where 

i  <  i’  and  k  <  k1 .  It  is  easy  to  show  that  the  covariance  matrix  %  is  identical 

2  P 

to  the  vector  having  K(K-l)/2  components  (o  -t)  *(Z  -Z  , )»  lSKi'SK,  where 

'  2  ?  , 

(Zx,  ....  ZR)  ~  NK(jO,IK),  a  =  E(G(Yi.))  =  %  and  t  =  E(G(Y1J)G(Yi,  .)) .  Q.E.D. 


Proof  of  Lemma  5.  The  inequality  D.  .  S  ot. -a.,  holds  if  and  only  if  at  least 

v*  /  i  i 

Z  of  the  differences  (Y^-a^)  -  (Y^,^,-a^,)  are  less  than  or  equal  to  zero  and 
n  n  2 

hence  1  Z  <t>(Y  -a,,  Y ,.,-a.,)<n  -2Z .  The  second  statement  of  this  lemma 

j=l  j'=l  1  1  J  1 


can  be  obtained  in  a  similar  manner. 


Q.E.D. 
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APPENDIX  B 


Players 


Rounding  first  base  times 
Methods 


Round  out  Narrow  Angle  Wide  Angle 


1 

5.40 

5.50 

5.55 

2 

5.85 

5.70 

5.75 

3 

5.20 

5.60 

5.50 

4 

5.55 

5.50 

5.40 

5 

5.90 

5.85 

5.70 

6 

5.45 

5.55 

5.60 

7 

5.40 

5.40 

5.35 

8 

5.45 

5.50 

5.35 

9 

5.25 

5.15 

5.00 

10 

5.85 

5.80 

5.70 

11 

5.25 

5.20 

5.10 

12 

5.65 

5.55 

5.45 

13 

5.60 

5.35 

5.45 

14 

5.05 

5.00 

4.95 

15 

5.50 

5.50 

5.40 

16 

5.45 

5.55 

5.50 

17 

5.55 

5.55 

5.35 

18 

5.45 

5.50 

5.55 

19 

5.50 

5.45 

5.25 

20 

5.65 

5.60 

5.40 

21 

5.70 

5.65 

5.55 

22 

6.30 

6.30 

6.25 
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